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ABSTRACT 

The contribution of Source Cosmic Rays (SCRs), accelerated and still con- 
fined in Supernova Remnants (SNRs), to the diffuse high energy 7-ray emission 
above 1 GeV from the Galactic disk is studied. The 7-rays produced by the SCRs 
have a much harder spectrum than those generated by the Galactic Cosmic Rays 
(GCRs). Extending a previous paper, a simple SNR population synthesis is con- 
sidered and the Inverse Compton emission from the SCR electrons is evaluated in 
greater detail. Then the combined spectrum of 7-ray emission from the Galactic 
Supernova Remnant population is calculated and this emission at low Galactic 
latitudes is compared with the diffuse 7-ray emission observed by the EGRET 
and ground based instruments. The average contribution of SCRs is comparable 
to the GCR contribution already at GeV energies, due to Supernovae of types 
II and lb exploding into the wind bubbles of quite massive progenitor stars, and 
becomes dominant at 7-ray energies above 100 GeV. At TeV energies the domi- 
nant contribution is from SCRs in SNRs that expand into a uniform interstellar 
medium. In fact, the sum of hadronic and Inverse Compton 7-rays would ex- 
ceed the limits given by the existing experimental data, unless the confinement 
time Tsni i-e. the time until which SNRs confine the main fraction of accelerated 
SCRs, is as small as T S n ~ 10 4 yr and the typical magnetic field strength in SNRs 
as large as 30 /iG. Both situations are however possible as a result of field ampli- 
fication through CR backreaction in the acceleration process. It is pointed out 
that accurate measurements of the low-latitude diffuse Galactic 7-ray spectrum 
at TeV-energies can serve as a unique consistency test for CR origin from the 
Supernova Remnant population as a whole. 
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1. Introduction 

The diffuse Galactic 7-ray emission was measured with the Energetic Gamma Ray 
Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory up to about 20 
GeV. To first approximation it can be described by a suitable model for the diffuse interstellar 
gas, the Galactic Cosmic Ray (GCR) distribution, and the diffuse photon fields (Hunter et 
al. f997a,b). The resulting GCR distribution in the Milky Way essentially corresponds to 
that measured in situ in the local neighborhood of the Solar System. This is also our starting 
point. 

However above about 1 GeV the observed average diffuse 7-ray intensity in the inner 
Galaxy, 300° < I < 60°, |6| < 10°, exceeds the model prediction significantly, and this excess 
appears to increase monotonically with 7-ray energy. Spatially, the diffuse 7-ray intensity 
is strongly increasing towards the Galactic midplane. 

In a previous paper Berezhko & Volk (2000a, hereafter referred to as Paper I), we have 
calculated the contribution of an unresolved distribution of CR sources to the diffuse 7-ray 
flux in the disk, assuming the GCR distribution within its large confinement volume to be 
the same as that measured in situ. The CR sources contribute through the 7-ray emission of 
the accelerated particles that are confined in their interior - the Source Cosmic Rays (SCRs). 
It was assumed that Supernova Remnants (SNRs) are the dominant sources of the GCRs. 
On this premise it was found that the SCRs give an important contribution to the diffuse 
high-energy 7-ray emission at low Galactic latitudes. Since the CR energy spectrum inside 
SNRs is much harder than on average in the Galaxy — the average spectrum being softened 
by rigidity-dependent escape from the Galaxy in the diffusion/convection region above the 
disk — the relative SNR contribution increases with energy and was found to become in fact 
dominant at 7-ray energies e 7 >100 GeV. This led to a substantial increase of the "diffuse" 
TeV 7-ray emission from the Galactic disk so as to constitute a significant observational 
background. This hard-spectrum background must also be taken into account in the search 
for spatially extended Galactic CR sources above the GeV region. 

Extending Paper I we take into account here that a significant fraction of the SNR 
population, the type II and lb Supernovae (SNe) with progenitor masses in excess of 15M , 
expand into a nonuniform circumstellar medium strongly modified by the wind from the 
progenitor star. In such SNRs the majority of CRs should be produced in a thin dense shell 
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of swept-up ISM. As a consequence of the high density of this shell the expected 7r°-decay 7- 
ray emission with energies e 7 < 100 GeV is shown to be an order of magnitude larger than for 
an identical explosion into a uniform ISM. This gives a natural explanation for the observed 
diffuse 7-ray excess for energies > 1 GeV. We also consider in detail the contribution of the 
high-energy SCR electrons to the high energy 7-ray background due to inverse Compton (IC) 
scattering on the background photon field. Since the energy spectrum of the SCR electrons 
is strongly influenced by synchrotron losses, the IC 7-ray emission depends on the effective 
magnetic field strength in the SNR interior. This field may in turn be amplified relative to 
the external field by the acceleration process itself (Bell & Lucek 2001). We take such a 
field amplification into account as an average property of SNRs relative to the contribution 
from the GCRs. As demonstrated below, this results in an order of magnitude increase of 
the Galactic 7-ray emission at TeV-energies. In a preliminary form these results have been 
presented in Berezhko & Volk (2003). 

We therefore conclude that a detailed measurement of the low-latitude diffuse Galactic 
7-ray spectrum within the energy interval from 1 to 10 4 GeV will constitute an important 
indirect test for the GCR origin from the population of Galactic SNRs as a whole. 

2. Gamma-ray luminosity of SNRs expanding into a uniform circumstellar 

medium 

The majority of the GCRs with charge number Z, at least up to kinetic energies e ~ 
10 15 Z eV, is presumably accelerated in SNRs. Individual examples, where this is the case, 
are the objects SN 1006 (Berezhko et al. 2002, 2003a) and Cas A (Berezhko et al. 2003b; 
Berezhko & Volk 2004). According to theory a significant part of the hydrodynamic SN 
explosion energy Esn ~ 10 51 erg is converted into CRs already in the early Sedov phase of 
the evolution, t ~ 10 3 yr, as a result of diffusive shock acceleration (Berezhko et al. 1996; 
Berezhko & Volk 1997, 2000b). Later on the total CR energy content E c and therefore the 
7-ray production slowly varies with time. We neglect this variation in the energetics in 
our considerations below, assuming a time-independent SCR energy content E c = O.IEsn, 
consistent with the GCR energy budget requirement. 

The total number of active SNRs, N S n = ^snTsn ~ 300 to 3000, is an increasing 
function of their assumed life time T S n = 10 4 to 10 5 yr, i.e. the time until which they can 
confine the main fraction of the accelerated particles. Here vsn ~ 1/30 yr is the Galactic SN 
rate. Consequently we conclude that the population of the oldest SNRs should dominate the 
total 7-ray luminosity of the ensemble of Galactic SNRs. Approximately we shall therefore 
consider only SNRs in the Sedov phase of their evolution in our estimate of the total 7-ray 
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luminosity. 

Gamma-ray producing CRs in the Galaxy are then the sum of two basically different 
populations. The first one consists of the ordinary GCRs and presumably occupies a large 
Galactic residence volume quasi-uniformly. This residence volume exceeds by far the volume 
of the gas disk which harbors the CR sources (Ptuskin et al. 1997; Berezinsky et al. 1990). 
The second CR population are the SCRs in the localized SNRs. 

During the initial, active period of SNR evolution at times t<T S N when the SN shock 
is relatively strong, the volume occupied by the accelerated CRs practically coincides with 
the shock volume. During later evolutionary stages the shock becomes weak and CRs begin 
to leave the SNR acceleration region. After some period of time the escaping SCRs become 
very well mixed with the ambient GCRs. Thus the controling factor is the shock strength. 

The 7-ray production by the GCRs is quite well studied (Hunter et al. 1997b; Mori 
1997). Therefore we shall concentrate on the relative contribution of the SCR population. 

It was shown in Paper I that the total 7-ray spectrum measured from an arbitrary 
Galactic disk volume is expected to be 

a p dF GCR 

^ = ^—[lA + R(e 1 )], (1) 

where dF GCR /de 1 is the 7r°-decay 7-ray spectrum due to GCRs and the additional fac- 
tor 0.4 is introduced to approximately take into account the contribution of GCR electron 
component to the diffuse 7-ray emission at GeV energies (Hunter et al. 1997b). 

For the ratio -R(e 7 ) = /Q^ of the 7-ray production rates due to SCRs and GCRs 
we have 

fl w = ^(ife)(n£v)°' 6 ( 1 + fl -»' < 2 > 

with £ = Ng CR /N GCR , where N GCR = lcm -3 and Ng CR are the gas number density in the 
Galactic disk and inside SNRs, respectively, R ep = Ql^ /Q pp is the ratio of the IC to 7r°-decay 
7-ray production rates due to the SCRs, T p is the proton SCR confinement time, and the 
difference between the spectral indices of the GCRs and the SCRs is assumed to be 0.6. 

In the case of SNRs expanding into a uniform ISM T p , which is equal to the confinement 
time T e of the SCR electron component if synchrotron losses are not more restrictive for 
electrons (see below), is given by the expression 

T p = mm{T SN , 10 3 ; (e / e max )~ 5 yr}, (3) 

where e max = 10 5 Z GeV is the maximal energy of CRs that can accelerated in SNRs, and 
e 7 = O.le which is roughly valid for the hadronic 7-ray production process. The decrease of 
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the proton confinement time T p with energy e is due to the diminishing ability of the SNR 
shock to produce high-energy CRs during the Sedov phase. In fact, the highest CR energy 
e and therefore the highest-energy 7-ray production slowly decrease with time as e oc t -1 / 5 
which implies the escape of the highest energy CRs from the SNR (Berezhko 1996; Berezhko 
et al. 1996; Berezhko & V61k 2000a). 

Since it was shown in Paper I that electron SCRs give a 7-ray contribution comparable 
with that of the nuclear SCRs, we consider the electron emission here in greater detail. In 
the relevant energy range the 7-ray luminosity due to the electron IC scattering on the 
background photons can be described in the Thompson limit (Longair 1981; Berezinsky et 
al. 1990). Since the electron spectrum in a SNR depends significantly on the SNR age t due 
to synchrotron losses, we first determine the 7-ray production rate dQlf (e 7 ) / dN$ N due to 
one SNR of age t: 

dQfi^) M de e dn e SCR (e e ) 

- °TcN ph — — , (4) 

where 

e e = m e c 2 y /, 3e 7 /(4e pA ) (5) 

is the energy of electrons which produce an IC photon with mean energy e 7 , Ot = 6.65 x 
10~ 25 cm 2 is the Thomson cross section, and e p h and N p h are the mean energy and number 
density of the background photons, respectively. Finally 

dn e SCR /dNl N = N e SCR /V g (6) 

denotes the spatial electron SCR number density averaged over the Galactic disk volume V g , 
where dN$ N = vsndt is the number of SNRs of age between t and t + dt which contribute 
to the IC emission at energy e 7 from CR electron sources. 

The source electron spectrum at energies e < where synchrotron losses are not 
important, can be represented in the form 

N e SCR (t,e) =K ep N SC R(e), (7) 

where K ep 10~ 2 is the electron to proton ratio and Nsciii^de is the total number of SCR 
protons per SNR in the energy interval de. Finally 

L25 . lo^wio^oy TeV 



t J \ B 

is the low energy limit of that part of the electron spectrum which is modified by the 
synchrotron losses, and B is the magnetic field strength inside the SNRs. For e > ei the 
electron spectrum is steeper due to synchrotron losses: 

N e SCR (t,e)=K ep N SCR (e)(e l /e). (9) 
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We also neglect the time dependence of the proton spectrum in SNRs, since in the later 
Sedov phase which is most important here, it is not significant, especially at high energies 
(Berezhko et al. 1996). The proton spectrum is taken in a power lower form Nscr oc e~ 7 
with 7 = 2.15, consistent with the requirements for the sources of the Galactic CRs. 

The maximum energy of SCR electrons is also restricted by their synchrotron losses 
during their acceleration at the SNR shock: 

V. \ ( B N - 1/2 



£ - = 24 UM^Jfel TCV ' <10) 

if e e max < e max ', here V s is the shock speed. At the beginning of the Sedov phase (t ~ 10 3 yr) 
the shock speed is about V s ~ 4 x 10 3 km/s and subsequently decreases as V s oc t~ 3 / 5 . 
Therefore the electron confinement time can be written in the form 

T e = min{T p , 10 3 (e e /6^J- 5 / 3 yr}. (11) 



" Te dO IC 



The total IC 7-ray production rate due to all existing SNRs can be determined by the 
expression 

/ dW-" SNdt (12) 

JO alV SN 

First we normalize Q^ c to the production rate of 7r°-decay 7-rays from inelastic CR - gas 
collisions, primarily p-p collisions, which may be written in the form (Drury et al. 1994) 

QP p (e) = Z^a pp cN g n(e), (13) 

where N g is the local gas number density, a pp is the inelastic p-p cross-section, Z 1 is the 
so-called spectrum-weighted moment of the inelastic cross-section, n(e) = N S cR(^)^SNTp/V g 
is the spatial number density of SCRs averaged over the Galactic volume V g , and c is the 
speed of light. Then, performing the integration over all possible SNR ages, we find the ratio 
of luminosities: 



Nph_\ f 4e 7 e pfe ^ 



(7"l)/2 

for 7-ray energies e 7 < e*, and 



/ NphT e \ f Ae^tph 



(7-l)/2 \ -1/2 



n^rssK^^j^) [n, + (18) 

for e 7 > e* where 




s 
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is the energy of 7-rays which are emitted by electrons with synchrotron loss time equal to 
T e . The value Z 7 = 0.113 was used, which corresponds to 7 = 2.15 (Drury et al. 1994). 

Due to the rather hard SCR electron spectrum the dominant contribution to the IC 
radiation comes from collisions with cosmic microwave background (CMB) photons which 
are characterized by N ph = 400 cm -3 and e ph = 6.7 x 10 -4 eV. The far infrared radiation 
(FIR) field with N ph « 20 cnT 3 and e ph « 0.01 eV contributes about 20% for e 7 > 100 GeV. 
In the case of the FIR background the Klein-Nishina cross-section must be used instead of 
the Thomson limit at TeV 7-ray energies. We take account of the Klein-Nishina decrease of 
the cross-section at high energies through multiplying a? by the correction factor 



which approximately describes the reduction effect (Blumenthal & Gould 1970). For e 7 = 
10 TeV, f KN m 0.1 for the FIR contribution. 

Bremsstrahlung 7-rays play no role for the average 7-ray background above GeV ener- 
gies, if K ep <§C 0.1 (Paper I). 

In Fig.l we present a calculated 7-ray spectrum based on the above expressions (1)— (17) 
with T SN = 10 5 yr, B = 10 fiG, e max = 10 5 GeV, K ep = 10~ 2 and N^ CR = Nf CR . A SCR 
power law index 7 = 2.15 and a 10 percent efficiency of SCR production in SNR, required 
to account for the GCRs, were used. 

One can see that the SCR contribution becomes dominant for e 7 >100 GeV where 
the expected 7-ray spectrum becomes extremely hard. At TeV energies the predicted flux 
exceeds the lowest HEGRA upper limit (Aharonian et al. 2001) almost by a factor of two. The 
SCR electron contribution increases with energy e 7 and exceeds the proton SCR contribution 
for energies e 7 ^300 GeV. At e 7 ~ 1 TeV electron SCRs contribute about 60% of the total 
7-ray flux. Therefore at these high energies the expected 7-ray emissivity of SNRs is only 
weakly dependent on the most relevant parameter, the maximum confinement time T$n, 
since according to expressions (2) and (15) we have for e 7 ~ 1 TeV 



taking into account that at this energy T p = T e = T SN . Even if T S n is as short as 3 x 10 4 yr, 
the expected SCR contribution at TeV energies is therefore only 30% lower and still exceeds 
the GCR contribution by more than an order of magnitude (see Fig.l). We emphasize that 
in the latter case the SCR contribution is entirely due to the electrons: they radiate about 
85% of the TeV 7-rays : This extreme case has to be considered as the lowest value for the 
SCR contribution to the Galactic 7-ray background radiation, since there are no other ways 




(17) 



dF* c /de^ oc 1 + 0.31n(TWl0 5 yr), 



(18) 
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to decrease it. This holds in particular also for the average gas number density Ng CR inside 
SNRs because the IC 7-ray production rate does not depend on it. Such a short confinement 
time T$n ~ 10 4 yr (Ptuskin & Zirakashvili 2003) could be realized in a scenario where a 
strong amplification of the magnetic field B in SNRs is balanced by nonlinear wave-wave 
interactions such that the wave magnetic power P{1) in the spatial scale / is dissipated in a 
(minimal) eddy turnover time ~ l/v(l) (Verma et al. 1996). Here v(l) = va[P(1) / 'B 2 ] 1 / 2 is 
the rms- value of the plasma mass velocity at scale I, and va denotes the Alfven velocity. 

We also point out that the restriction of the electron confinement time T e due to their 
synchrotron losses becomes significant for 7-ray energies e 7 > 1 TeV and leads to the steep- 
ening of the 7-ray spectrum (see Fig.l). 

According to Fig.l, for energies e 7 = 0.4 to 4 TeV the expected differential 7-ray energy 
spectrum from the central part of the Galaxy has the form 

dF^/de-f = A(e 7 /1 TeV)"" (19) 

with power law index a « 2.25 and amplitude A = (6 to 8) x 10 -9 (cm 2 s sr TeV) -1 , 
depending on the confinement time T SN . 

As one can see from Fig.l, the 7-ray flux expected for T S n = 10 5 yr at e 7 = 1 TeV 
exceeds the HEGRA upper limit (Aharonian et al. 2001) and the preliminary flux measured 
by the Milagro detector (Fleysher 2003). We emphasize here that we have multiplied the 
Milagro flux by a factor of 3 since we present the expected 7-ray flux from the region |6| < 2° 
whereas the Milagro data correspond to \b\ < 5°. Even for T$n = 3 x 10 4 yr the expected 
flux exceeds the Milagro data. 

The only physical parameter which strongly influences the IC 7-ray production rate is 
the SNR magnetic field B: according to Eqs. (15)— (16) oc B~ 2 . A value of B that is 
significantly higher than B = 10 /xG can be attributed to field amplification at the shock 
front due to the strong wave production by the acceleration of CRs far into the nonlinear 
regime (Bell k Lucek 2001). 

In Fig. 2 we present the same calculations as in Fig.l but with a mean SNR magnetic 
field value B = 30 fiG. Compared with the previous case the IC 7-ray emission is decreased 
by an order of magnitude. Therefore the contribution of the hadronic SCRs to the 7-ray flux 
becomes much more significant: at TeV-energies and for T$n — 10 5 yr the 7r°-decay 7-rays 
exceed the IC 7-rays by a factor of 2.5 , whereas for T S n = 3 x 10 4 yr the IC 7-rays still 
exceed the 7r°-decay 7-rays by a factor of 1.3. One can see that in this case the expected 
7-ray flux is below the HEGRA and Tibet upper limits already for a SCR confinement time 
T S n = 10 5 yr. The expected 7-ray spectrum at energies e 7 = 0.1 to 10 TeV is characterized 
by a = 2.4 and A = (2.4 to 4) x 10 -9 (cm 2 s sr TeV) -1 , which also considerably exceeds 
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the GCR contribution. Note that since the maximum energy of SCRs obeys e max oc B, it 
is larger by a factor of three compared with the previous case. Taking into account the 
uncertainty of the Milagro flux one can conclude that this measurement does not contradict 
our prediction for the case of the amplified magnetic field. The lower value of the confinement 
time Tsn ~ 10 4 yr might be preferable in this sense. 

3. Wind SNRs of type II and lb 

The above consideration corresponds to the assumption that all existing Galactic SNRs 
expand into a uniform circumstellar medium with a density of N^ CR = 1 cm -3 . However, for 
stellar masses in excess of about 15 M the progenitors of type II and lb SNe can strongly 
modify their environment through their ionizing radiation and very energetic winds (Weaver 
et al. 1977; Chevalier & Liang 1989). Disregarding large-scale turbulent mixing processes 
in the sequel, the stellar wind from such a star creates a low-density bubble beyond the 
termination shock out to a radius R s h, bounded by a dense swept-up shell of interstellar 
material. The radius of the wind bubble can be several tens of parsecs. The very hot dilute 
bubble contains a small amount of mass, typically a few solar masses or even less, whereas 
the thin swept-up shell ultimately contains thousands of solar masses. This shell will have a 
gradual inner boundary where localized mixing of bubble and shell material proceeds. Since 
the typical SN ejecta mass is about M e j pa 10M Q for massive progenitors, only a relatively 
small fraction of the SN explosion energy E$n is transformed into internal gas and CR 
energy during the SN shock propagation through the bubble and the main part of SN energy 
is deposited in the shell. 

There are no physical reasons for which the CR production process in the shell matter is 
substantially different from that in the uniform ISM case, apart from the different distribution 
of ion injection across the SNR shock surface due to the magnetic field geometry (Volk et al. 
2003). The major difference which plays a role for the 7r°-decay 7-rays is that CRs accelerated 
during the SN shock propagation through the shell are confined in a medium which is more 
than a factor of ten denser than the average ISM in the Galactic disk, of density N^ CR . This 
may lead to an appreciable increase of the overall 7-ray production. 

In addition, the SNR shock evolves differently than in the case of a uniform circumstellar 
medium. The entire active period of SNR evolution, when the SN shock is strong enough and 
effectively produces CRs, takes place within the shell of thickness L which is much smaller 
than its radius, L <C R s h- This implies that during the entire active evolution the SNR shock 
size R s remains nearly constant. The time dependence of the shock speed V s which is very 
important for the SCR confinement in SNRs can then be estimated as follows. 
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The majority of the progenitor stars with an intense wind are main-sequence stars 
with initial masses Mj > 15M Q (Abbott 1982). In the mean, during their evolution in the 
surrounding uniform ISM of gas number density po — m p N^ CR , they create a bubble of size 
(Weaver et al. 1977; Chevalier & Liang 1989) 

R sh = 0.76(0.5MV; 2 4/po) 1/5 , (20) 

where M is the mass-loss rate of the presupernova star, V w is the wind speed, and t w is the 
duration of the wind period. 

In order to determine the SN shock dynamics inside the shell we model the gas number 
density distribution in the bubble and in the shell in the form 

N g = N b + (r/R sh ) 3 ^N sh , (21) 

where N sh = aN^ 011 is the peak number density in the shell, N b is the gas number density 
inside the bubble, typically very small compared with the shell density, and a = N s h/N is 
the shell compression ratio. 

The mass of the bubble 

M b = (A7TR 3 sh /3)m p N b (22) 
is about 3 to 5M , whereas the shell mass 

M sh = 4nN sh m p I drr 2 (r / R sh ) 3 ^ = (4irR 3 sh /3)N° CR m p (23) 
Jo 

is a few thousands of solar masses. Therefore, during SNR shock propagation through the 
bubble, only a small fraction of its energy is given to gas of stellar origin. The main part of 
the explosion energy is deposited in the shell. 

As in the case of a uniform ISM the evolution of the SN shock can be represented as a 
sequence of two stages. 

The first stage is the free expansion phase which continues up to the time, when the 
shock sweeps up a mass equal to the ejecta mass M e j. The second phase is analogous to the 
Sedov phase in the uniform ISM. The only significant difference compared with the uniform 
ISM case is that at the end of the free expansion phase, when a large part of the SN energy 
has been given to the swept up gas, the shock volume is much larger. 

When the swept-up mass becomes large compared to the sum of bubble mass and M ej , 
the SN shock propagation in the shell medium with density N g oc r 3 ^ -1 ) can be shown to 
approach the following self-similar adiabatic solution 

R 3 s V s 2 N g (R s ) = const, (24) 
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which corresponds to an expansion law R s = R (t/to) u with v = 2/ (3a + 2). The shock 
decelerates much more rapidly in this case, V s oc t~ 3a ^ 3cT+2 \ than in the case of a uniform 
ISM, where V s oc t~ 3 / 5 . Due to the fact that 3d > 1 we have approximately V s oc t _1 . During 
this period of time the shock size and the amount of the swept-up mass are connected by 
the relation 

R s = R sh {M/M sh f'^. (25) 

Substituting in this expression the SN ejecta mass M ej we get the value of the SN shock 
size which corresponds to the end of free expansion phase 

R s (t ) =R = R sh (M sh / M ej )~ 1 ^ 3a . (26) 

At this epoch the shock speed is roughly the mean ejecta speed 

V a mV = ^2E SN /M ej . (27) 

Therefore the corresponding time scale is t = Rq/V . 

As it follows from the expansion law and from eq.(24), the SNR shock reaches the edge 
of the shell R s (tf) = R s h at time 

t f = t (M sh /M ej f° +2 ^, (28) 

which for a 3> 1 gives 

*/ « t ^M sh /M ej . (29) 

At about this stage the SNR shock has lost much of its speed and will come into pressure 
equilibrium with its environment, i.e. tf ~ T S n- 

For the typical values M = 6 x l(T 9 M /yr, V w = 2500 km/s, t w = 10 7 yr, which 
correspond to a progenitor of mass M,- L = 15M Q (Abbott 1982; Chevalier & Liang 1989) and 
n gcr = 1 cm -3 ; we have Rgh K go pc anc i M sh = 3x 1O 3 M . For E SN = 2 x 10 51 erg 

and M e j = 5M Q this gives to — 5 x 10 3 yr and tf ~ 10 5 yr. The ambient interstellar gas 
may also be photoionized which implies a large enough external pressure to stop bubble 
expansion at a lower radius (Chevalier & Liang 1989). Therefore we also consider below the 
case R s h = 17 pc which corresponds to tf = 3 x 10 4 yr. 

According to theory, the expanding SNR shock produces a power law CR spectrum up 
to a maximum energy (Berezhko 1996; Berezhko et al. 1996; Berezhko & Volk 1997) 



e m oc R S V S , 



(30) 
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which is determined by the radius R s and speed V s of the shock. The CRs with the highest 
energy e max are produced at the very beginning of the Sedov phase t ~ t when the product 
R S V S has its maximum. Subsequently, the product R S V S decreases with time approximately 
as t^ 1 and the SNR shock produces CRs with progressively lower cutoff energy e m (t) < 
tmax = Cm(^o)- During that phase those CRs that were previously produced with energies 
e m < e < e max now leave the remnant without a significant influence of the SNR shock. 
Therefore the CR confinement time has to be taken in the form 

T p = mm{l, (t /t f )(e max /e)}T SN , (31) 

where according to the above consideration T$n — tf- This relation is analogous to eq.(3) 
in the case of a uniform ISM and leads to an increase of the spectral index of the 7-ray 
emission by one unit. 

In order to estimate the highest energy of accelerated CRs in this case we use the 
relation e max oc RqVqB. Compared with the case of a uniform ISM, where Rq = 4 pc for 
N g 3CR = 1 cm -3 , the value of Rq is about ten times larger but Vq is the same. The magnetic 
field B(to) on the other hand is much lower than in the case of a SN la. Indeed, as follows 
from eqs. (25) and (20), the gas number density at the beginning of the Sedov phase is as 
small as N g (t ) = N sh M ej /M sh w N^ CR /300. The magnetic field in the region of the shell is 
presumably adiabatically compressed interstellar field B . If we approximate the directions 
of Bq as being isotropic, then its magnitude scales with the gas density like B oc Ng^ 3 , 
e.g. Chevalier (1974). This gives a field strength B(t ) 0.02£? . Taking into account all 
the factors considered, we conclude that the maximum CR energy in the case of wind-SNe 
with dense shells is roughly a factor of ten smaller than in the case of SNe expanding into 
a uniform ISM. Assuming that the field amplification is proportional to y/N~ g , cf. Bell & 
Lucek (2001), this amplification does not play a role here since N g (to) is so low. 

In contrast to our previous study (Berezhko & Volk 2000b), where we considered the CR 
acceleration by SN shocks expanding into a modified bubble, we consider here the opposite 
extreme of bubble structure. It is assumed that the magnetic field suppresses the mass 
and heat transport between the dense shell and the hot bubble (Chevalier & Liang 1989). 
Therefore the bubble mass is so small that a significant number of CRs is produced by the 
SN shock only when it enters the dense shell. Even though it is not clear at the moment 
which of these two concepts is physically more correct, we consider here the possibility of an 
unmodified bubble, because of its great importance for 7-ray production. 

We note that due to the energy dependence of the confinement time T p (e) the SCR 
contribution to the 7-ray flux undergoes a spectral break at a 7-ray energy which corresponds 
to the SCR energy 

^break 

(t /tf)e max . (32) 
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Protons with energies e < e br eak survive up to the final active SNR epoch t — tf and therefore 
provide the maximum contribution to the 7-ray SNR emissivity at e^ eak = 0.1e break , whereas 
for e 7 > e b J eak the 7-ray emissivity decreases more steeply to go to zero at e 7 = 0.1e max . 

The gas number density N g SCR in expression (2) is a function of SNR age t for wind-SNe, 
because the gas density N g (R s ) seen by the SN shock changes during the evolution of the 
SNR. It then follows from eq. (20) that the gas density behind the SN shock, where most of 
the SCRs are located, can be represented as N 2 (t) = <JsNN g (R s (t)) = <JsN(t/tf) 2 N s h, where 
(Tsn is the SN shock compression ratio. Therefore the gas density seen by SCRs of energy e 
is 

N s g CR = a SN (T p /T SN ) 2 N sh . (33) 

Since T p oc e^ 1 for e 7 > e^ eak , the 7-ray spectrum produced in the shell is very steep at 
these high energies, and in eq. (2) we have i2(e 7 ) oc e 7 2 ' 4 - 

We assume that the progenitors of type II SNe are stars more massive than Mj > 8M , 
and that only those with Mj > 15M have strong winds which produce extended bubbles, 
e.g. Abbott (1982). This is also true for the progenitor population of type lb SNe. According 
to Giisten & Mezger (1983) the initial mass function has the form dN/dMi oc M^ a with 
a = 1.6,2.4,3.24,3.62 for the mass intervals Mj < 1M , 1M < Mj < 1OM , 1OM < 
Mj < 5OM , 5OM < Mj, respectively. Therefore the ratio of the number of stars whose 
initial mass Mj exceeds 15 M to those with Mj > 8M is 0.23. Since about 85% of all SN 
explosions in the Galaxy are type II and type lb SNe (Tammann et al. 1994), and since in 
addition the explosion energies E$n of type II and type lb SNe with Mj > 15M are about 
two times larger than on average (Chevalier 1977; Hamuy 2003), we conclude that roughly 
a fraction 5 = 0.3 of the total Galactic SN energy release is from wind-SNe with progenitors 
of masses Mj > 15M which produce extended bubbles. 

Therefore the expected 7-ray SNR luminosity should be weighted between the SNR 
populations expanding into a uniform ISM and into bubbles as follows: 

R = (l- 5)i? 7 + SR n , (34) 

where Rj corresponds to the ratio of 7-ray production rates due to SCRs and GCRs if all 
SNe were exploding into a uniform ISM, and Rn is the corresponding ratio for explosions 
into bubbles. 

Assuming that wind-SNe have rarefied bubbles with dense shells the spectrum of 7-rays 
calculated with the above ratio R, a = 10, a SN = 5 and 5 = 0.3 is presented in Figs.l and 2 
by the full lines, which correspond to two different assumptions about the value of the SCR 
confinement time T S n- It is seen that the wind-SNe dominate at energies e 7 < 100 GeV and 
fit the EGRET data for T S n = 3 x 10 4 yr fairly well, whereas for the larger confinement 
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times T S n = 10 5 yr their contribution becomes so large, that their calculated flux exceeds 
the EGRET data by about a factor of two. 

We emphasize that the contribution from the population of wind-SNe is mainly due 
to the nuclear SCR component. The break in the 7-ray spectra presented in Fig.l, 2 is 
at e^ eak = 30 GeV for T SN = 10 5 yr and at e^ eak = 100 GeV for T SN = 3 x 10 4 yr. For 
larger energies the proton confinement time and the gas density become lower, T p oc e _1 , 
Ng oc e~ 2 , which leads to a corresponding decrease of the wind-SN contribution to the 
7-ray spectrum. This contribution becomes insignificant for e 7 ~ 1 TeV. 

4. Discussion 

Our above consideration is based upon a picture in which the population of SNRs, 
considered here as the main GCR sources, are distributed across the Galactic disk similar 
to the Interstellar gas. In this case the calculated ratio R{e 1 ) of the SCR to the GCR 
contribution does not depend upon the observing direction in the disk. Since the number 
of these sources N$n — ^snTsn is very limited, their mean number within the field of 
view ~ few degrees of a stereoscopic system of imaging atmospheric Cherenkov telescopes 
for 7-ray energies >100 GeV is so low, that one should expect large fluctuations of the 
actual value of i?(e 7 ), especially for lines-of-sight other than those directed towards the inner 
Galaxy. 

The actual SNR distribution within the Galactic disk most probably is not uniform. Ac- 
cording to Case & Bhattacharya (1998) the SN explosion rate as a function of Galactocentric 
radius r has a peak at r 5 kpc and drops exponentially with a scale length of ~ 7 kpc. 
Within 5 < r < 20 kpc this agrees fairly well with the radial distributions of Supernovae in 
a sample of 36 external galaxies of Hubble type Sb-Sbc, statistically considered equivalent 
to our Galaxy (Dragicevich et al. 1999). From our results, for 7-ray energies larger than 10 
GeV, the ratio R, corrected for the nonuniform SNR distribution, is then given by the ratio 
R c = i](r)R, where the parameter r](r) represents the ratio of the actual SNR number along 
the given line-of-sight to that SNR number which would corresponds to their uniform distri- 
bution. For lower 7-ray energies, 100 < e 7 < 10 4 MeV, where the truly diffuse emission from 
the GCRs dominates, it is known that the 7-ray emissivity gradient is significantly shallower 
than the gradient of rj(r), as shown by the EGRET instrument on board of the CGRO satel- 
lite (Strong & Mattox 1996) in basic agreement with the results from the COS-B satellite 
(Strong et al. 1988). This much shallower, indeed truly diffuse Galactocentric 7-ray gradient 
can be understood as a nonlinear propagation effect from the disk into the Galactic Wind 
where the latter is driven by the GCRs themselves (Breitschwerdt et al. 2002). In contrast, 
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the radial galactocentric variation of the direct radiation from the sources is not smoothed 
by any propagation effect. Therefore its amplitude should spatially vary like the distribution 
of the sources, integrated along the line-of-sight, and this variation should become directly 
measurable in the > 100 GeV range. Given several such line-of-sight integrals, one can infer 
the radial galactocentric 7-ray emissivity gradient. The excellent angular resolution of atmo- 
spheric Cherenkov telescopes and the high sensitivity of instruments like the CANGAROO 
III and H.E.S.S. arrays in the Southern Hemisphere, which at TeV energies have an order of 
magnitude higher sensitivity compared with previous generation instruments like HEGRA 
(Aharonian et al. 2001), should allow this measurement in particular in the inner Galaxy. 

The assumed average value B = 30 /iG of the magnetic field inside SNRs is considerably 
smaller than the value 120 fiG inferred for e.g. SN 1006 in the very early Sedov phase 
(Berezhko et al. 2002, 2003a). However the amplification must be a monotonically increasing 
function of the shock strength. In particular, Bell & Lucek (2001) argue for a dependence 
B cc V s . If so, then the effective amplified magnetic field in SNRs varies from B ~ 100 /xG 
at the beginning of the Sedov phase to the typical interstellar value B pa 10 fiG at the end 
of its active period t = Tsn- Since the value of the magnetic field in SNRs has a dominant 
influence on the IC 7-ray spectrum, and since all evolutionary phases of a SNR contribute 
to this spectrum in a roughly equal way, the adopted value B = 30 /iG appears to be quite 
realistic, averaging over the entire active period. 

One should also note that the SNR magnetic field influences the maximum energy of 
SCRs. Since the most energetic SCRs are created at the beginning of the Sedov phase, 
their actual maximum energy e max oc B(to) is expected to be larger than considered here, 
because the amplified magnetic field B(t ) in this phase is by a factor of several larger than 
on average. Compared with the Fig. 2 this effect will produce more a extended and smooth 
high energy tail of the 7-ray spectrum up to the energy e 7 10 5 GeV. 

The actual process of SN shock interaction with the dense shell could be much more 
complicated compared with the ideal picture considered here. If due to some physical fac- 
tors (inhomogeneities in the surrounding ISM, instabilities ...) the shell is always strongly 
distorted, then one would expect that the SN shock will penetrate through the shell more 
rapidly. This will lead to a decrease of the propagation time tf and/or to the decrease of the 
effective gas number density Ng CR , that in turn will decrease the amount of 7-rays produced 
in such types of SNRs. In other words, the calculated 7-ray spectrum at e 7 ^ 100 GeV (see 
Fig.l, 2), which is mainly due to the wind SNe, has to be considered as an upper limit, which 
corresponds to the case of stable unmodified bubbles. 
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5. Summary 

Our considerations demonstrate that the SCRs inevitably make a strong contribution 
to the "diffuse" 7-ray flux from the Galactic disk at all energies above a few GeV, if the 
population of SNRs is the main source of the GCRs. According to our estimates, the SCR 
contribution dominates at energies greater than 100 GeV due to its substantially harder 
spectrum. 

There are two physical parameters which influence the expected 7-ray emission from 
SNRs significantly: the CR confinement time and the mean magnetic field strength B inside 
the SNRs. For a conventional value B = 10 /xG the expected 7-ray flux from SNRs exceeds 
the HEGRA upper limit considerably if the SCR confinement time is as large as T SN = 10 5 yr. 
This contradiction can be resolved either if we suggest an appreciably higher postshock 
magnetic field B>30 /xG or if the SCR confinement time is as small as T$n ~ 10 4 yr, or a 
combination of both effects. These possibilities can be attributed to field amplification by the 
SCRs themselves. In fact, nonlinear field amplification may also lead to a substantial decrease 
of the SCR confinement time: according to Ptuskin & Zirakashvili (2003) maximal turbulent 
Alfven wave damping with its corresponding increase of CR mobility could make the SCR 
confinement time as small as T S n ~ 10 4 yr. Under these circumstances the most realistic 
7-ray background spectrum is represented by the thick solid line in Fig. 2. One can see that 
even in the case of B = 30 /xG and T SN ~ 3 x 10 4 yr the SNR contribution at TeV energies, 
with roughly numbers of IC and 7r°-decay 7-rays ; still exceeds the GCR contribution almost 
by an order of magnitude. Note that the preliminary 7-ray flux at e 7 = 1 TeV, measured 
by the Milagro detector (Fleysher 2003), confirms our earlier prediction (Berezhko & Volk 
2003). 

At lower energies e 7 ;$ 100 GeV the 7-ray emission from SNRs is dominated by the wind- 
SNe with initial progenitor mass M« > 15M Q which expand into the bubble created by the 
progenitor's wind, under the assumption that the bubble is not modified by global heat and 
mass transport. The main fraction of the CRs is produced in this case when the SN shock 
propagates through the thin dense shell at the edge of the bubble. Due to the significantly 
higher mass density of the shell material compared with the mean gas density in the Galactic 
disk 7r°-decay 7-rays dominate at 7-ray energies e 7 < 100 GeV despite the fact that only about 
20% of all SNRs belong to this class of objects. As shown in Figs.l and 2, the discrepancy 
between the observed "diffuse" intensity and standard model predictions at energies above 
a few GeV can be attributed to the SCR contribution. This requires unmodified bubbles for 
the wind-SNe. Since we cannot prove this assumption, the corresponding explanation of the 
lower energy 7-ray excess is not a definitive conclusion, but rather a plausible suggestion 
which must await more detailed studies of wind bubble morphologies. 
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We conclude that a detailed measurement of the low-latitude "diffuse" Galactic 7-ray 
spectrum within the energy interval from 1 to 10 4 GeV will allow a strong consistency check 
for the predominant origin of the Galactic Cosmic Rays from the Galactic population of 
SNRs as a whole. In the energy range above a few 100 GeV our predictions are quite robust 
resulting in a hard spectrum whose Galactocentric variation should correspond to that of the 
observed SNR distribution. If the preliminary Milagro data, which agree with our prediction, 
are confirmed by more precise measurements, then such measurements will give an indirect 
confirmation of SNRs as the main sources of GCRs. The measurements will also give a new 
tool to study the spatial distribution of SNRs in the Galaxy. 
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Fig. 1. — The average diffuse 7-ray spectrum of the low latitude inner Galaxy (38° < I < 43°, 
\b\ < 2°). The thick and the thin dashed lines represent the case of a uniform ISM with the 
two different SNR life times T S n = 10 5 yr and T S n = 3 x 10 4 yr, respectively; the full lines 
correspond to the case when 30% of the SN energy release is due to SNRs which expand 
into the wind bubble of a massive progenitor star. EGRET data (Hunter et al. 1997b), 
preliminary Milagro data at e 7 = 1 TeV (Fleysher 2003), the Wipple upper limits (Reynolds 
et al. 1993; LeBohec et al. 2000), the HEGRA upper limits (Aharonian et al. 2001), and the 
Tibet upper limits (Amenomori et al. 2003) at e 7 = 3, 10 TeV are shown. The Milagro data 
have been multiplied by a factor of 3 (see text). 
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